Clustering Properties of Dynamical Dark Energy Models 
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We provide a generic but physically clear discussion of the clustering properties of dark energy 
models. We explicitly show that in quintessence-type models the dark energy fluctuations, on scales 
smaller than the Hubble radius, are of the order of the perturbations to the Newtonian gravitational 
potential, hence necessarily small on cosmological scales. Moreover, comparable fluctuations are 
associated with different gauge choices. We also demonstrate that the often used homogeneous 
approximation is unrealistic, and that the so-called dark energy mutation is a trivial artifact of an 
effective, single fluid description. Finally, we discuss the particular case where the dark energy fluid 
is coupled to dark matter. 



I. INTRODUCTION 

Observations in the last decade have gradually been 
providing evidence suggesting that not only most of the 
universe's energy is in a dark, unknown form, but indeed 
the dominant component of this dark sector violates the 
strong energy condition — the latter being required to ex- 
plain the recent acceleration of the Universe. From a 
purely phenomcnological point of view, the simplest can- 
didate for this task is Einstein's cosmological constant 
(leading to the so-called ACDM or concordance model), 
and indeed this is in decent agreement with most of 
the existing data. However, from a fundamental physics 
point of view, the vacuum energy density suggested by 
observations is many orders of magnitude smaller than 
the most optimistic expectations, so that models where 
dark energy is dynamical (for example, being due to a 
scalar field) are arguably the more likely ones. 

While there is, in some sense, a single constant dark 
energy model, the number of possible dynamical dark 
energy models is clearly infinite. Moreover, the standard 
observational techniques in use for probing dark energy 
are not ideal, in the sense that they do not measure di- 
rectly the theoretically relevant quantities, and they al- 
low for several important degeneracies to remain when 
comparing theory and observation. Both of these mean 
that even in an era of precision cosmology there is ample 
room for building phenomenological models of increasing 
complexity, whose differences when it comes to cosmo- 
logical observables can be quite small, and hence are not 
easy to distinguish even with good data. 

One may legitimately ask if these are merely 'epicycles' 
which arc doomed to be replaced by something entirely 
different and much simpler. While a fairly strong case 
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can legitimately be made for a positive answer to this 
question, what should replace them still remains to be 
determined. What one can certainly say at this stage is 
that increasingly elaborate mathematics should not be 
used as an excuse to obscure simple physics. In other 
words, if one builds more elaborate models on the basis 
of pre-existing paradigms, then the broad general fea- 
tures of the original paradigm will still be present and 
can seldom be avoided, even if there are differences in 
the specific details. 



The goal of the present paper (which is a follow-up and 
in the same spirit of [H) is to present a physically clear 
discussion of the clustering properties of dynamical dark 
energy models, specifically those where the dark energy 
is due to a classical scalar field with a perfect fluid form. 
We shall be relating our discussion to other recent lit- 
erature, contrasting the different approaches in order to 
clarify crucial issues. We will provide an explicit simple 
proof that in quintessence-type models the dark energy 
fluctuations, on scales smaller than the Hubble radius, 
are of the order of the perturbations to the Newtonian 
gravitational potential, and hence necessarily small on 
cosmological scales. We show that the often used homo- 
geneous approximation is unrealistic and does not pro- 
vide useful information, at least on sub- horizon scales. 
Moreover, we illustrate two physically simple points that 
are often obscured by unnecessary mathematics: the fact 
that fluctuations comparable to the Newtonian gravita- 
tional potential can be obtained by a mere gauge trans- 
formation, and that the so-called dark energy mutation 
is a trivial artifact of an effective, single fluid descrip- 
tion. Finally, we discuss the particular case where there 
is a coupling between the dark energy and dark matter 
fluids. Throughout this paper we shall work in natural 
units and a metric signature (— , +, +, +). 
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II. SCALAR FIELDS 



We shall be interested in the action 



S = J d A x y^g£{X, . 



(1) 



where C is the Lagrangian for a real scalar field (f> and 



(2) 



is the kinetic term. We shall generically consider the 
case of a classical scalar field <fi governed by the an ar- 
bitrary Lagrangian of the form £(X,<f>). We can write 
the energy-momentum tensor for this model in a perfect 
fluid form 
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On the other hand, the scalar field equation of state is 
defined by 



P 

w = - 



C 



f(X) - V{<t>) 



2XC,x - C 2Xf, x -f + V^) 



(7) 



Observe that by carefully choosing V(4>) we may inde- 
pendently specify w and c 2 . This simple point will be 
important for some of what follows. Finally, we can also 
differentiate Eq. ([7]) and substitute Eq. ([6]) to find 



w.x = {c 2 s - w) 



e .x 



(8) 



which shows that w = c s is sufficient to ensure that the 
equation of state has no dependence on X. It is not, 
however, necessary as it also occurs whenever the energy 
density itself has no such dependence (e.x = 0); the cos- 
mological constant is a simple example. 



by means of the following identifications 
V u < 



In Eq. ([3]), u M is the 4- velocity field describing the mo- 
tion of the fluid (for timelike V^), while e and p are 
its proper energy density and pressure, respectively. Ob- 
serve that from this it trivially follows that if p — p(X), 
then e = e(X). Unfortunately it is not always possible 
to invert e(X) and obtain X(e), but when it is the fluid 
has an explicit isentropic equation of state p = pis). 

In the special case where C(X, <p) = f(X) — V(<f>) (with 
/ and V being arbitrary functions of X and 0, respec- 
tively), the mass of the scalar field is defined as 



8 2 V 



(5) 



For a viable dark energy model one typically requires a 
mass smaller than the Hubble parameter H , m < H , at 
least for a canonical scalar field with f(X) = X. On 
small (effectively sub- horizon, < m" 1 ) scales, pressure 
and density fluctuations are related through Sp = c?5p 
where the sound speed is given in linear theory by [2| 



c 2 = 



p,x 
e ,x 



C 



C 



X 



2XC 
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(6) 



If f(X) — X we have a canonical scalar field which has 
a constant sound speed of unity. However, we have the 
freedom to choose the function / in order to obtain very 
different sound speeds. An example of an algebraically 
simple but physically interesting class of functions is 
f(X) cx A™, which yields c 2 = l/(2n - 1). In partic- 
ular, when n = 1, the scalar corresponds to a massless 
scalar field, n — 2, to background radiation, and so on. 
In the limit n — > oo, the scalar can be interpreted as dust 
(in other words, pressureless non-relativistic matter). 



III. HOMOGENEOUS APPROXIMATION 

Many aspects of the evolution of dark energy pertur- 
bations have been studied in great detail, particularly in 
the context of standard quintessence scenarios. These 
are characterized by two main components, dark matter 
(DM) and dark energy (DE), and the latter is usually 
modeled as a canonical scalar field with f(X) = X. As 
shown in the previous section, dark energy in these mod- 
els is characterized by c 2 — 1 (and w ~ —1)- However, 
in many of these studies, only the homogeneous case has 
been investigated [3J, U, l5( . For our present purposes, the 
key feature of this approximation is that pressure gradi- 
ents are completely absent. 

Although the homogeneous approximation greatly sim- 
plifies the analysis, it is not difficult to see, even in the ab- 
sence of any detailed comparisons, that the correspond- 
ing results are unrealistic. Indeed, such an approxima- 
tion can only be considered reliable when uniform per- 
turbations on scales larger than the Hubble radius are 
considered, which is clearly not the situation of greater 
physical interest. Still, before proceeding to our main 
analysis, we will briefly review the homogeneous case (in 
which dark matter and dark energy collapse together), 
so as to provide a useful comparison point. 

Let us assume that the energy-momentum tensors of 
DE and DM are separately conserved, in other words, 
that there is no coupling between them (we will discuss 
the coupled case later in this note). If we consider the 
evolution of a region of physical volume V, we can write 



so that 



cfeDE + (1 + H^DE^J- = 



dV 

«£dm + £DM-y- — U , 



<5de = (1 + w)S-dm , 



(9) 
(10) 

(11) 
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where w = pde/ede is the dark energy equation of state, 
5 = (e — e)/e is the contrast and e denotes the average 
density. 

Simple as this is, it is sufficient to highlight an im- 
portant consequence of this approximation. For generic 
values of the equation of state parameter w (viz., values 
not incredibly close to w = —1), large fluctuations in the 
dark matter component will inevitably be accompanied 
by significant fluctuations in the dark energy component. 
As we will see in the next section, this is in fact not the 
case for realistic sub- horizon perturbations. 



IV. LINEAR FLUCTUATIONS 

An order-of-magnitude estimate of the perturbations 
in the dark energy component can be obtained by con- 
sidering linear theory. We will therefore consider a stan- 
dard dark matter fluid with wdm = and a dark energy 
fluid with and arbitrary equation of state w and sound 
speed c 2 . Unless otherwise stated we shall be working in 
the synchronous gauge, which is comoving with the dark 
matter, that is u BU = 0. 

The linear evolution of scalar perturbations is de- 
scribed by 



where we have taken the sound speed c 2 = 1 and also 
defined M = Se^mL 3 - This shows that any fluctuations 
in the dark energy component are of the order of the per- 
turbations to the Newtonian gravitational potential and 
thus necessarily very small on cosmological scales. This 
result agrees with those of @, 0] , but has been obtained 
in a physically clearer way. 

Indeed, with hindsight this result is hardly surprising. 
A simple way to understand it is to observe that in the 
standard quintessence scenario, in order to have w ~ 
— 1, one requires |X/V| <C 1. It then follows that the 
amplitude of the fluctuations on sub-horizon scales must 
necessarily be small in this type of models. For this not 
to be the case, the kinetic term X would wave to be 
substantial and we would no longer have a dark energy 
fluid. The only way to get around this restriction is to 
consider other classes of dynamical dark energy models. 

However, this is not yet the full story. It must be 
emphasized that fluctuations of this same order are also 
associated with different gauge choices. For example, by 
considering a gauge transformation from a local inertial 
frame comoving with the dark matter fluid (which we will 
denote as frame I) to the local inertial cosmological frame 
(denoted by frame II, with a vanishing CMB dipole) , one 
obtains 



$t>m + T~18dm — - 7~i 2 [^dm<5dm + (1 



3c 2 )JWde] = 
(12) 



6de + 3H(c 2 - w)S DE + (1 + w){9 - 5dm) = (13) 



9 + H(l-3c 2 s )9 



-V 2 6 



1 + w 



DE 







(14) 



where a dot denotes an rj conformal time derivative 
(with drj = dt/a and TL = aH = a/a), il; = e\je c 
(e c = 3H 2 /8irG being the critical density), and we have 
also defined 



9 = V ■ ude = oKe), 



(15) 



On scales much smaller than the Hubble radius (that is 
kTi S> 1), and assuming for simplicity that the dark mat- 
ter fluctuations are the only source for the dark energy 
fluctuations, we can write 



£de,ii — £de,i 
£de,i 



wv 2 D 



and this can be approximately written 



A 



(l + w)v 2 D 



(1 + w) 



GM 



(19) 



(20) 



where Vjj is the velocity of the dark matter fluid with 
respect to frame II and 7d = (1 — v 2 D )~ x l 2 . In the last 
step we have made the reasonable assumption that the 
main contribution for the dipole is directly related to 
the local Newtonian gravitational potential induced by a 
mass perturbation M with length scale L. Perturbations 
of comparable magnitude may be similarly obtained by 
considering other possible gauge choices. 



DARK ENERGY MUTATION 



1 + w 



-V 2 <5 



DE 



-H 2 O dm c>dm = 47rGa 2 <5e D M • (16) 



If we now consider a density perturbation of comoving 
wave number k one immediately finds that 



c 2 fc 2 £pE 
1 + w 



4irGa 2 5e 



DM ■ 



(17) 



Equivalently, we can write this as a function of the char- 
acteristic scale of the perturbation L ~ a/k, 



6 



DE 



GM 



1 + w 



(18) 



The analysis of the previous section shows that if the 
sound speed of the dark energy fluid is significant, then 
the dark energy perturbations must necessarily be small. 
However, for very small values (c 2 ~ 0) the above state- 
ment is no longer true, so this case warrants a separate 
treatment. In this section we shall therefore consider a 
model with a null sound speed c 2 = 0, by requiring the 
pressure p to be constant, and study the behavior of a 
region of physical volume V. 

Energy-momentum conservation trivially implies 



d[(e + p)V} = 0, 



(21) 
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and consequently (e + p) oc V 1 . 
can be written in the form 

1 



The equation of state 



(22) 



where C is a constant. This simple calculation shows that 
w will depend on the physical volume V of the region 
under consideration and, consequently, on the density 
perturbations. This effect is just a coarse-graining issue, 
in the sense commonly discussed in condensed matter 
systems. We see that in low-density regions, where V is 
very large, w ~ — 1, but in collapsed regions the value of 
w can be much smaller (in modulus) and even approach 
zero: indeed, in the limit V — > 0, we have w — ► 0. 

The most noticeable consequence emerging from this 
analysis is that if density perturbations are present, then 
it necessarily follows that w cannot be a constant. This 
effect has been pointed out in Q and dubbed 'dark energy 
mutation'. However, there is nothing surprising about 
it, and indeed its physical explanation is rather prosaic. 
This model turns out to be the simplest example of uni- 
fied dark energy models studied in [J, [l(| , and it has 
been shown to be totally equivalent (to any order) to 
ACDM. In fact, there is no way to distinguish the single 
fluid interpretation made above from the standard inter- 
pretation with canonical components, dark matter and a 
cosmological constant, having respectively u>dm = and 
wde = —1. Dark energy mutation, in this case, is just 
an artifact of our single fluid description. 



VI. COUPLED MODELS 

In the so-called concordance model, a range of obser- 
vational data is used to postulate the existence of two 
dark fluids (DE and DM) for which so far there is no 
direct experimental backing. In the context of GR, the 
most common attitude is to model DE and DM as two 
minimally coupled fluids. The direct opposite to this is 
to treat them as different manifestations of a single fluid 
(UDE models). An intermediate approach, on the other 
hand, is to view them as coupled fluids. In this case, 
however, if the coupling is very strong, we naturally ex- 
pect the distinction of DE and DM as two different fluids 
to become somewhat blurred. In other words, we expect 
that, to a certain extent, strongly coupled fluids will be- 
have as if a single fluid. As far as we are aware today, 
this bridge between strongly coupled models and UDE 
has not been significantly explored. 

Following the recent literature [E 0, ED, 03, 0, 
[lU , we shall assume that the two dark fluids are coupled 
through the Lagrangian 



C = X - V{(f) + h((j>)£ 



DM 



(23) 



(cj) in this context is normally called a 'chameleon' field 
[111, [IE IS EH] •) Note that according to the discussion in 
Sect. UH the dark matter component may be described 



by a scalar field ip with a Lagrangian 

£dm oc Y n , 
in the limit of large n, where 



(24) 



(25) 



Indeed, for a fixed value of n, <p> has an equation of state 
parameter wum = (2n — I) -1 and, therefore, pdm = 
wdm£dm becomes negligible for large n. If follows that 
we can rewrite (|2"B1 as 



(26) 



L = X - V{<f>) + . 9 (</>)£DM 



where g(4>) = wdmM'W i s a rescaled coupling constant. 
It is easy to check (by varying the action in relation to tp) 
that the dark matter component evolves independently 
from the chameleon field. On the other hand, the evolu- 
tion of 4> is given by 



U<f> = 



dV cS 



where 



VeS = V((f>) - g{4>)ST>M ■ 



(27) 



(28) 



and therefore is affected by how dark matter evolves. 
(Note here that although V c s ~ V, dV c s/d<t> can be very 
different from dV/d(f>.) As for the energy-momentum ten- 
sor associated with (1231) , it is a simple matter to show (by 
varying the action in relation to g^) that 

T^fa tp) = V M 0V„0 +{X- V(<t>)) giu , + 

+ h(^)[g^Y n +nY n - 1 V^pV l ,tp] .(29) 

Obviously, this energy- momentum tensor does not, in 
general, describe a perfect fluid. However, in the so- 
called adiabatic regime (described, in detail, in [l4l [l5|). 
it is assumed that the gradients of 4> are both negligible 
in T^ v and the equation of motion (|2T[) . Thus, in this 
regime, reduces to 



T M „ ~ (hY n - V)g^ + nhY n ~^^V vV 



(30) 



which can be immediately written in a perfect fluid form, 
if we make the following definitions 



^ V2Y 

PcS = -V e f[ 



£cff = ^DM 

-V. 



v. 



The effective equation of state is then 



Poff 

e c ff 



-v(4>) 



V(4>) + he 



DM 



(31) 



(32) 



Is it an isentropic fluid, though? Yes. Since the adiabatic 
regime is also characterized by the condition dV Q g /dip = 
0, the value of cj> is univocally related to e c s. Hence, 
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Peff only depends on the value of e e g and therefore the 
fluid is isentropic (although, in general, we won't have an 
explicit p e f{ = Pes(£es) equation of state). 
Now, the value of 



where ' denotes d/d<fi. It follows that 



d 2 V cS 



(33) 



determines the length scales for which the adiabatic 
regime is valid. Specifically, this is the case for large 
scale perturbations with L , while for scales much 

smaller that this, the approximation is no longer valid. 
We thus conclude that above a certain scale, sufficiently 
coupled models behave as a single isentropic fluid but 
not below. Why is this relevant? It is relevant because 
the background evolution in UDE models described by a 
single isentropic fluid is expected to be strongly affected 
by non-linear effects which severely complicate the anal- 
ysis of this type of models. On the other hand, it is 
still unclear if the differentiated behaviour above or be- 
low L ~ "m-efi may help to relax the averaging problem 
(Iol | that affects UDE models. For the sake of argument, 
suppose the majority of the non-linear clustering occurs 
for scales smaller that m~g ; since now they are confined 
to a non-isentropic part of the fluid, it is possible that 
they may not greatly affect the large scale evolution of 
the universe. If, on the other hand, significant clustering 
does extend beyond this scale, then non-linearities will 
still be a major problem in strongly coupled models (at 
least for models not sufficiently close to a ACDM model) . 
Finally, note that at recent epochs, 



TO) | » |<?(#?dm|, 

and, in the adiabatic regime 



(34) 



(35) 



v{4>) 



< 



5(0) 



G'{4>) 



G(0) 



(36) 



Given the stringent astrophysical and laboratory con- 
straints on variations of Newton's constant G, the above 
relation imposes strong constraints on the shape of the 
potential in the region explored by the field in recent 
times. 



VII. CONCLUSIONS 



In this paper we studied in some detail the clustering 
properties of dynamical dark energy models. We have 
shown that in standard quintessence models, contrary to 
recent claims, dark energy fluctuations on sub-horizon 
scales are of the order of the perturbations to the Newto- 
nian gravitational potential and, consequently, very small 
on cosmological scales. We have also pointed out that the 
homogeneous approximation is not an adequate frame- 
work to study the evolution of dark energy fluctuations 
on small scales. 

We have also explored the extent to which coupled 
models can be interpreted as unified models. Our anal- 
ysis makes it clear that the non-linear instability that 
plagues unified dark energy models [l[ should also apply 
to these so-called chameleon models. A more detailed 
analysis of this scenario is warranted, but we leave it for 
future work. 
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